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Using a Green's function approach, we investigate band 
structure effects in the BEEM current distribution in recip- 
rocal space. In the elastic limit, this formalism provides a 
"parameter free" solution of the BEEM problem. At low 
temperatures, and for thin metallic layers, the elastic approx- 
imation is enough to explain the experimental I(V) curves 
at low voltages. At higher voltages inelastic effects are ap- 
proximately taken into account by introducing an effective 
RPA-electron lifetime, much in similarity with LEED theory. 
For thick films, however, additional damping mechanisms are 
required to obtain agreement with experiment. 
PACS numbers: 61.16.Ch, 72.10.Bg, 73.20.At 

Ballistic Electron Emission Microscopy (BEEM), and 
its spectroscopic counterpart (BEES) were originally 
designed as techniques extending the power of Scan- 
ning Tunneling Microscopy (STM) to buried interfaces, 
particularly of metal-semiconductor systems. The stan- 
dard model describes BEEM as a convolution of three 
steps Q: tunneling from the tip (1), propagation in the 
metallic layer (2) and transmission through the metal- 
semiconductor interface (3). This model clearly suggests 
the important potential of BEEM to focus in any of these 
steps separately. However, it is unnecessary to stress that 
such a deconvolution process may only be safely per- 
formed applying a sufficiently elaborated theory, which 
should use as few adjustable parameters as possible. In 
the past, the lack of such a precise method to analyze the 
experiment has prompted several intense discussions: (i) 
whether fcy is conserved or not at the interface Q, (ii) 
the origin of the observed nanometric resolution and its 
relation to the tunneling injection [Q, (iii) the similar re- 
sults obtained on Au/Si(lll) and Au/Si(100) interfaces, 
despite of their different projected conduction-band min- 
ima 1^ , (iv) how ballistic are the electrons in BEEM af- 
ter all? 0, etc. This list of intensively debated questions 
in the literature is probably an indication of the limita- 
tions associated with the standard approach based on E- 
space Monte-Carlo simulations, where processes crucial 
from a physical point of view are simply parametrized 
to give agreement with experiment. In particular, in all 
these Monte-Carlo calculations, the energetic spectrum 



and the momentum distribution of the injected electrons 
are taken from conventional planar tunneling theory, us- 
ing a free electron approach. This assumption is proba- 
bly the origin of the major limitation to a first-principles 
analysis of the BEEM current, as the propagation of the 
electrons in the metal film is strongly dependent on the 
metal band structure and can depart significantly from a 
free electron behaviour [||. Accordingly, the aim of this 
letter is to present a microscopic formalism that incorpo- 
rates those band structure effects and yields the appropri- 
ate angular momentum distribution that, as shown below 
for the case of gold films, drastically departs from the 
narrow forward cone assumed in E-space Monte-Carlo 
simulations. We will show how most of the previous in- 
terpretations of BEEM data for Au/Si interfaces need to 
be modified when using the right k-space currents. 




FIG. 1. (a) k-space current distribution for Au(lll) inside 
the coherence region (5th atomic layer) {V = 1 eV, = 0.1 
eV); (b) same as (a), but outside the coherence region (30th 
atomic layer). Dark regions correspond to higher intensities. 
The Si BZ (small) and the Au BZ (large) are shown together 
with the ellipses where the Si conduction band minima project 
(notice that the outer ellipses appear after the corresponding 
remapping) . 

We introduce a full quantum-mechanical description 
of the BEEM problem based on a Keldysh Green's func- 
tion method written in a Linear Combination of Atomic 
Orbitals (LCAO) basis. Our analysis is based on the 
following three-step scheme: We provide an accurate de- 
scription of the initial tunneling injection (1), and the 
subsequent propagation of electrons through the metalic 
layer (2). Passing over the Schottky barrier (3) is taken 
into account applying energy and fcy conservation, and 
matching states at the two-dimensional interface. In this 
paper the foregoing scheme is applied to the case study 



of a (111) oriented gold metallic layer deposited on a 
(111) silicon substrate; applications to other metals (e.g. 
CoSi2) and other semiconductors are in progress. Per- 
fect unrelaxed surfaces and bulk-like ideal geometries are 
assumed in our analysis, but it is seen from the nature 
of our results that a relatively small amount of disorder 
(e.g., confined to 3-4 layers close to the interface) would 
not fundamentally change our conclusions. 

A Green's function formalism presents the important 
advantage of being free of any adjustable parameter in 
the strictly elastic limit, where we only add an arbitrarily 
small positive imaginary part to the energy (77), necessary 
to ensure attenuation of the wave at infinity. Moreover, 
inelastic effects associated with the electron-electron in- 
teraction can be added incorporating a complex energy 
dependent self-energy, 7]{E). We shall view the self- 
energy as a single parameter to be adjusted to the experi- 
ment, representing an effective inelastic electron-electron 
mean free path producing attenuation: Xatt ~ -^^^ ■ This 
method has been succesfully adopted to different fields, 
like Low-Energy Electron Diffraction. 

In an LCAO basis, we write the Hamiltonian as: 



H — Ht + Hs + Ht 



(1) 



where Ht — ^afia + Tapc]^Cfj defines the tip (Greek 
subindices), Hs = Yl^i^i + S^y^^i designates the 
metal substrate (Latin subindices), and Hj = ^ TajC^^Cj 
describes the coupling between the tip and the metal sur- 
face in terms of a hopping matrix, Taj, expressed as a 
function of the different atomic orbitals in the tip and 
the surface by using a tight-binding formalism (n^, 
cjj, and Cq, are number, creation and destruction opera- 
tors defined in the usual way). 

Since the system under investigation is out of equilib- 
rium, a convenient way to compute the current between 
two sites i and j in real space is given by Keldysh's tech- 
nique P|: 



J^,= j Tr{f,,{Gt--G%-)]dE 



(2) 



The matrices G^" are non-equilibrium Keldysh Green's 
functions that can be calculated in terms of the standard 
retarded and advanced Green's functions 10|. We no- 
tice that this formalism allows us to compute on the same 
footing the tunneling current between the tip and the 
sample and the current propagating in the metal (steps 
1 and 2). To this point, all our expressions are exact, 
and the main task is to determine how to compute the 
retarded and advanced Green functions, and which ap- 
proximations are introduced there. 

Previously |^ , we have analyzed the electron propaga- 
tion in real space, using a semiclassical approximation for 
these Green functions, and have found important focus- 
ing effects in gold films. Now, we concentrate on calculat- 
ing the full quantum-mechanical current distribution in 



reciprocal space, using a formalism based on rcnormaliza- 
tion group techniques . This momentum distribution 
will allow us to obtain the spectral I(V) characteristics. 
In particular, the current between two layers a and h 
inside the metal, at a given energy E and fc|| , can be 
expressed as [| 12 : 



•MA) I 



(3) 



where g^^^ '{E,k\\) is the retarded (advanced) Green's 
function for the unperturbed metal linking the layer h 
and the surface layer, 1, Tab{k\\) is a hopping matrix con- 
necting layers a and &, and poo{E) is the density of states 
on the last atom of the tip (0), considered for simplicity 
to be the only tip active atom for tunneling. The trace 
denotes a summation over the orbitals forming the chosen 
basis. 

Step three of our four-step scheme involves comput- 
ing the transmission coefficient for the two-dimensional 
interface. Applying a surface Green's function match- 
ing formalism jl^ ] in the neighbourhood of the M point, 
we obtain a transmission coefficient T{E, fcy) that can be 
used in k-space to give the injected current in the semi- 
conductor: 

pEp+eV /• 

HV) = dE dk\\J,_i^,{E,k\^) X T(^,A:||) 

JEp+eVo JlstB.Z. 

(4) 

where c refers to the metal layer at the interface, and Vb 
is the Schottky barrier height (assumed to be 0.86 eV); 
note that the transmission coefficient is zero outside the 
ellipsoids allowed by energy conservation (see Figure 1). 
The integral inside the first Brillouin zone is performed 
summing over a dense grid of special points |l4[ . 

In previous publications we have discussed how the 
propagation of electrons in the gold periodic lattice re- 
sults in focused beams and narrow Kossel-like lines in real 
space, with a 3-fold symmetry associated to the (111) di- 
rection of an fee crystal j|jl^,|l^. These lines have typ- 
ical widths of around 3-4 atomic distances, explaining 
the nanometric resolution of the BEEM technique even 
in deeply buried interfaces. These results also show how 
the Bloch wave is formed after propagation by more than 
four or five layers, forbidding the propagation of electrons 
in gap directions over longer distances. We believe that 
our results are convincing enough to answer a question 
nowadays found in the literature related to Monte-Carlo 
simulations: is it realistic to assume that electrons can 
propagate 20 or 30 A as free particles along the forbidden 
Au(lll) directions? We conclude that this is an unphys- 
ical scenario because of the strong deflection exerted by 
the lattice on electrons traveling in these directions. 

However, in this work we shall focus on our results in 
reciprocal space and their influence on the I(V) curves. 



An important feature observed in reciprocal space is a 
change in the symmetry of the k-space current distribu- 
tions when going from thin to thick layers. The expected 
symmetry for a quantum-mechanical calculation is re- 
lated to the projected density of states ||T^. It is six- fold 
in (111) fee planes, because of the equal contribution of 
+k and —k states. This is indeed the case for an arbi- 
trarily small imaginary part (77) added to the energy, but 
as commented above, rj can be interpreted in terms of 
a complex self-energy arising from inelastic events defin- 
ing a coherence region of the order of Xatt ■ Beyond that 
region inelastic processes become important, and inten- 
sities rather than amplitudes add to give the final wave- 
field. This takes us from a quantum mechanical picture 
(six-fold) to a semi-classical one (three- fold), as can be 
seen comparing Figure la, inside the coherence region, 
to Figure lb, where the current distribution is computed 
in a layer outside that area. The three-fold symmetry is 
progressively built up as a function of metal thickness, 
and can be understood in terms of our previous analysis 
the symmetry of the wavefield in the semiclassical 
limit is related to the Fermi surface, reflecting the three- 
fold symmetry of the crystal. Therefore, this is a new 
example of how a quantum system, under the influence 
of friction, becomes gradually classical by a decoherence 
process [T^ ]. In addition, it is seen how the current in 
k-space deviates for these thick layers from a simple den- 
sity of states calculation |l^, concentrating around the 
directions predicted by the semiclassical analysis (Figure 
lb) [pp^. The difference observed in reciprocal space 
between the quantum and semiclassical regime does not 
significantly affect 
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FIG. 2. Theoretical I(V) curves for Au/Si(lll), d=75 
A (experimental values -solid circles- from ref. 17): ballis- 
tic -dotted-, RPA approximation for \att{E) (see Eq.5) with 
a" — 260 A(eV)'^ -dashed line- and with Xatt modified -solid 
line- are shown. In the inset these two Xatt{E) are displayed. 

the beams in real space (where the symmetry must 
always be three-fold), but could in principle affect the 
I(V) current injected through the projected ellipses into 



the semiconductor. However, because of the gradual 
crossover seen from one regime to the other, we do not 
expect dramatic effects, unless one could experimentally 
break the time-reversal symmetry suddenly (e.g., by ap- 
plication of a magnetic field), or could selectively block 
the current injected in some of the six equivalent ellip- 
soids. In those sudden jump between a semiclas- 
sical regime and a quantum one should be observable. 
It should be noted that band structure effects results in 
k-space current with enough fcy to have the electrons in- 
jected into the outer conduction-band minima of Si (the 
central one is forbidden because of the gap in that direc- 
tion), and explains the long standing puzzle of why the 
threshold on Au/Si(lll) and Au/Si(100) is nearly the 
same: our calculations show how the similar results ob- 
tained for both interfaces are related to these nontrivial 
distributions in k-space, after the appropriate folding of 
the gold Brillouin zone inside the silicon one is performed 

Next, we compute theoretical I(V) curves from Eq. (2). 
A quantitative comparison with BEES experiments [ pT| , 
will then allow us to discuss also the electronic mean 
free paths. First of all, we try the hypothesis of ballis- 
tic electrons. On intuitive grounds this should suffice for 
low temperature, low voltages, and very thin layers. In 
Fig. 2 we compare experimental results for Au/Si(lll) 
at T=77 K, d=75 A with a pure ballistic calcula- 
tion (77 very small and injection at first attempt). It 
is clear from these results that, without using any ad- 
justable parameter, the onset is reasonably explained by 
a purely ballistic theory that uses the right current distri- 
butions in k-space. Therefore, we are able to give a rea- 
sonable explanation of the experiment for voltages near 
the threshold, but it is also noticed in Fig. (2) that data 
beyond V=1.2 eV can only be consistently interpreted 
by assuming an attenuated wave. To introduce an atten- 
uation mechanism, people have considered three major 
sources of damping: electron-electron, electron-phonon 
and electron-defect interaction. As the electron-phonon 
contributions are greatly reduced at 77 K we first con- 
sider a \att{E) dominated by the electron-electron inter- 
action. Within an RPA approximation for a free electron 
gas with a density representing gold (r^ — 3.01), we ob- 
tain: 



E/Ef 
{E - Epf 



(5) 



with A*^ = 260 A (eV)^. Results considering multi- 
ple reflections |jl^ between the surface and the interface 
through a specular model are presented in Fig. (2), where 
an excellent agreement is seen again up to 1.2 eV. Beyond 
that voltage, a reduction in the attenuation length by 
about 20% on average is required to bring experimental 
and theoretical intensities close together. The resulting 
Xatt[E) is displayed in the inset of Fig. (2). The reduc- 
tion with respect to the first-principles RPA approxima- 



tion might be understood as representing either band- 
structure or impurity effects in the effective eiectron- 
eiectron interaction. It is remarkable the good agreement 
obtained for low voltages, where the Xatt changes quickly 
with energy following an RPA-like behaviour. This is at 
variance with E-space Monte-Carlo simulations whereby 
a smoother dependence of Xatt with energy was found 
pTf . Our results suggest, however, that for thin films and 
low voltages, the main source of damping is the electron- 
electron interaction that is well described within a RPA 
approach. 
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FIG. 3. Theoretical I(V) curves for Au/Si(lll), d=300 A 
(experimental values -solid circles- from ref. 17): RPA with 
A° = 175 A(eV)^ -dashed line-, and with Xatt = 125 A for all 
E -solid line-. 

However, a different example of BEES-data, where a 
pure ballistic theory is not sufficient even near the thresh- 
old, is afforded by the case of thick layers (see Fig. 3). 
In this case we notice that if we use a RPA-like energy 
dependence for Xatt, we find both a discrepancy in magni- 
tude and a different voltage dependence for I(V) (as seen 
in the different slopes). If we choose a different A° in the 
RPA expression to get the right magnitude, we still would 
observe a serious discrepancy with the experiment (e.g., 
see Fig. (3) where A" has been reduced to 175 A (eV)^). 
Because all the other elements in the theory that might 
be responsible for the discrepancy {J{E) and T{E) in 
formula (4)) are calculated from first principles, we take 
this as a serious indication of a different dependence of 
Xatt{E) with E. A possible physical origin for this effect 
is the likely presence of defects (e.g. vacancies) The 
natural choice for this scenario is an energy-independent 
attenuation length in the Green function. With this as- 
sumption we obtain an excellent agreement with the ex- 
periment (T = 77 K, d = 300 A 0) for Xatt = 125 A, as 
seen by the solid line in Fig. 3. This value is in reasonable 
agreement with attenuation lengths derived by different 
groups in films of similar thickness |17 1^ and suggests a 
different behaviour of Xatt{E) between thick (300 A) and 
thin (75 A) films. 

In conclusion, we have introduced a Green's function 



formalism that in the ballistic limit is an ab initio ap- 
proach to BEEM. The particular k-space current dis- 
tributions determined by band structure effects are the 
main result of our analysis and crucial for a quantitative 
comparison with experimental BEEM data. Inelastic ef- 
fects have also been approximately included by use of 
an imaginary self-energy. This single quantity is fitted 
to the experiments to explain a number of spectroscopic 
data on the Au/Si interface. 
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